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Abstract — High orders of perturbation theory can be calculated by the Lipatov method, whereby they are 
determined by saddle -point configurations, or instantons, of the corresponding functional integrals. For most 
field theories, the Lipatov asymptotic form has the functional form ca N Y{N + b) (N is the order of perturbation 
theory) and the relative corrections to it are series in powers of UN. It is shown that this series diverges facto- 
rially and its high-order coefficients can be calculated using a procedure similar to the Lipatov one: the Kth 
expansion coefficient has the form const[ln(5 1 /5 )]~^r(/T + (r { - r )/2), where S Q and S { are the values of the 
action for the first and second instantons of this particular field theory, and r Q and r 1 are the corresponding num- 
ber of zeroth-order modes; the instantons satisfy the same equation as in the Lipatov method and are assumed 
to be renumbered in order of their increasing action. This result is universal and is valid in any field theory for 
which the Lipatov asymptotic form is as specified above. © 2000 MAIK "Nauka/Interperiodica ". 



1. INTRODUCTION 

Lipatov proposed a general method of calculating 
high orders of perturbation theory whereby these are 
determined by saddle-point configurations, or instan- 
tons, of the corresponding functional integrals [1]. On 
its appearance, the Lipatov method stimulated major 
discussion (see the collection of articles [2]) but was 
subsequently cast into doubt because of the possible 
existence of additional renormalon contributions [3]. In 
a recent work [4], the present author put forward a 
detailed discussion of the existing argumentation in 
support of renormalons and showed that this is untena- 
ble in the broad philosophical sense and in the mathe- 
matical sense: this clears any obstacles from applying 
the Lipatov method to a wide range of problems in the- 
oretical physics. 

The Lipatov method can be used to study any quan- 
tities [5] but the starting point is that it can be applied 
to functional integrals having the form 



AM) 



Kg) = jDcpcp (i) ..y 

xexp(- S {y}-gS int {y}), 



(1) 



where cp (1) , cp (M) is a certain sample q>,- , cp,- 

from the set of integration variables cp, contained within 
the symbol Dep. The expansion coefficients I N of the 
integral (1) in terms of the coupling constant g are 
determined by the Cauchy integral 



2ni g N+1 ' 



(2) 



in which the saddle -point method can be used for large 
N. The functional form of the Lipatov asymptote is 
given by 



I N = ca Y{N + b), N ^oo, 



(3) 



and the relative corrections to it have the form of a reg- 
ular expansion in terms of 1/7V: 



I N = caT(N + b) 
l + 4l + ^ + ...+4£ + 



(4) 



N 



AT 



N* 



Calculation of the corrections to the asymptotic form 
provides important information on the expansion coef- 
ficients and is an alternative to the direct diagrammatic 
calculations of the lower orders: for instance, instead of 
calculating the fourth or fifth orders [6, 7], it is more 
economical to calculate A! or A 2 . So far, the first correc- 
tions to the asymptotic form have only been calculated 
in cp 4 theory [8] and in a few quantum-mechanical prob- 
lems [9, 10]. 

In the present paper, we study the behavior of the 
coefficients A K for large K. This topic has not been stud- 
ied theoretically and the only available data has been 
obtained by numerical methods: for a perturbation the- 
ory series in the problem of an anharmonic oscillator, 
Bender and Wu [9] determined the first ten coeffi- 
cients A K : 

A y = -1.3194444, A 3 = -7.0142876, 
A 2 = -1.9385609, A 4 = -40.118943, 
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Fig. 1. (a) Saddle points and integration contour in integral 
(11). (b) In calculations of the asymptotic form of A K the 
contour must be deformed since the point t = corresponds 
to a singularity not a saddle point. 



A 5 = -305.5223, 
A 6 = -2808.09, 



A 8 = -3.65 x 10 , 



A 9 = -4.4 x 10 , 



(5) 



A 7 = -2.995 x 10 , A w = -lx 10 . 

The rapid growth of these coefficients indicates that the 
series in (4) diverges. 

Another example which can easily be studied is the 
zero-dimensional limit of cp 4 theory. In this case, the 
functional integral in fact reduces to a single one 



1(g) = p(pcp M exp(-(p 2 - <? (p 4 ), 



(6) 



and its expansion coefficients can be calculated in the 
explicit form: 



In = 



rl ,, M + 3V7,, m+ 1 



27271 



r(N + 1) 



(-4f . (7) 



By isolating the asymptotic form for N — - °° and 
expressing the result in the form (4), we obtain the fol- 
lowing expression for the coefficients A K for K — - °o 
(see Appendix) 



9 ( 1 -i- KiM \ 
A K = Re ZU+ ^ + / r(^), 



(8) 



whose functional form is similar to the Lipatov asymp- 
totic form (3) but with the complex parameters a and c. 

In the present study, we shall show that factorial 
divergence of the series in (4) also occurs in the general 
case and a universal result [see formula (47)] valid for 
any field theory using the Lipatov asymptotic form (3) 
can be obtained for the asymptotic form A K . 

2. SIMPLE EXAMPLE 
AND QUALITATIVE PICTURE 

The qualitative aspects involved in the calculation of 
the asymptotic form of A K can be conveniently demon- 



strated by calculating the corrections to the Stirling for- 
mula: 



T(N + 1) = j2nNe~ N N N 
xIi+4i + ^ + ... + 4£ + ... 



(9) 



N 



The result is well known for the logarithmic form of 
expression (9): in this case it is possible to find a gen- 
eral term of the series known as the Stirling series [11]. 
By calculating the exponential function of the Stirling 
series using factorial series algebra [5], we can easily 
find the asymptotic form of A K : 



A K = -Re 



2r(K) 

.jr+i' 



K 



(10) 



(271/) 

We shall subsequently show how this result is obtained. 

Using the definition of a gamma function and mak- 
ing the substitutions x — ► Nx and t = lnx, we have 

nN+\) = \dxx N e x = Ne- N N N 



x jfi?xexp{-/V[x- 1 - lnx]} (11) 

o 

= Ne~ N N N jdte'exp{-N[e' -l-t]}. 



For large N, the saddle-point condition has the form 
e* - 1 = so that there is a set of saddle points t s = 2%is, 
s = 0, ±1, +2, . . ., lying on an imaginary axis (Fig. la). 
The integration contour in (1 1) passes through the sad- 
dle point t = and satisfies all the conditions for valid- 
ity of the saddle-point method [12]. Thus, its deforma- 
tion is not required and the other saddle points can be 
neglected. Calculating the integral in the saddle-point 
approximation yields the Stirling's formula. 

Formally isolating the asymptotic form, we identi- 
cally set 



r(N+l) = j2%Ne~ N N N F(l/N) 
and making the substitution 
e = UN, 

we have for the function F(e) introduced by us 



/27te J 



dte exp< — 



1 



(12) 



(13) 



(14) 
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Expanding (14) as a series in terms of e gives the 
required coefficients A K which are calculated by anal- 
ogy with (2): 



de 1 




(15) 



-Kkie. 



This is an exact expression which for large K can be 
calculated using the saddle -point approximation. The 
saddle-point conditions yield the set of solutions 

t s = 2nis, e s = ~, s = integer, (16) 

so that in the complex plane t the saddle points are for- 
mally the same as those in the calculations of the lead- 
ing asymptotic term. However, for the integrand at the 
5th saddle point we can easily obtain the estimate 



~exv{-K-Klne s } 



tfK\, 



(17) 



from which it is clear that the solution with s = does 
not in fact correspond to a saddle point but to a singu- 
larity. 1 Hence, the integration contour over t cannot 
pass through the point t = but must be deformed and 
pass through one of the neighboring saddle points, 2ni 
or -2ni (Fig. lb) which, because of (17), gives the 
required asymptotic form of A K ~ (2n)~ K K\ [see (10)]. 2 

A similar situation is encountered in the general 
case. When the coefficients A K are calculated in the sad- 
dle-point approximation, the instanton equation is the 
same as that used to calculate the Lipatov asymptotic 
form. However, using the same solution as in the last 
case yields a singularity rather than a saddle point 
(because of the additional integration over e). There is 
thus a need to consider other solutions of the instanton 
equation which can be numbered in order of increasing 
action corresponding to them. If the Lipatov asymp- 
totic form is determined by the first instanton, having 
the smallest action, the principal contribution to the 
asymptotic form of A K is made by the second instanton. 

3. GENERAL CASE 

Calculations of the Lipatov asymptotic form (3) are 
fairly cumbersome if the aim is to find all its parameters 
a, b, and c. However, if the analysis is confined to deter- 



1 The solutions (16) are written assuming 6^0 which does not hold 
for s = 0. A similar observation must be made with reference to 
formula (38) below. 

2 The integration contour over e in (15) is conveniently drawn 
slightly to the right of the imaginary axis, enveloping the left half- 
plane over an infinitely distant contour; in this case for line < 
the integration contour over t is shifted upward and passes 
through the point 271/ whereas for line > it is shifted downward 
and passes through the point -2ni. 



mining the parameters a and b, it is possible to have 
simple structural calculations which reduce to a formal 
expansion near the saddle point and isolate the depen- 
dence on N. We shall demonstrate these calculations for 
the case of cp 4 theory; however, we do not need the 
explicit form of the action and we shall only use its 
characteristic properties of homogeneity 

5 (Xcp) = X 2 S Q (q>), S int a<?) = ^S int (<p). (18) 

Similar properties of homogeneity are encountered in 
other field theories and, with slight modifications, the 
scheme put forward subsequently also holds for the 
general case. 

According to (1) and (2), the expansion coefficients 
are given by 



/„., - M§w..v" 

c 

xexp(- S {<p}-gS int {<p}-N\ng). 
We introduce the new variable 

and set 

5{cp} = S {q>} + S int {$}. 



(19) 



(20) 



(21) 



In terms of the new variable (]), the saddle-point condi- 
tions have the form 



S'lA} = o, gc = 



N ' 



(22) 



and expanding the expression in the exponential func- 
tion (19) as far as quadratic terms in 8(|) = (]) - (p c and 
og = g - g c gives 

N(H, 5"{(b r }8(b) N , 2 
-N-N\ng c - ^ y ' V7 ^ - A ;(5gf. (23) 
2 S{<M 2g c 

We use a symbolic notation, denoting the first and sec- 
ond functional derivatives by single and double primes 
and taking these to be a vector and a linear operator, 
respectively; the variables of integration cp, contained 
within Dcp are taken to be components of the vector cp. 
Bearing in mind that because of (20) 



5c> = Jg~ c [S§ + ^-y c }, 5cp = cp-cp e , 



(24) 



and shifting the origin Sep, we have 

T -N -N+l-M/2 f dt f _ . (1) , 

xexp(-i(5(p, S"{<t> c }5q>) + f/), 
where we have set hg = ig c t. 



(25) 
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We make the linear substitution 8cp — ► 5 Sep with 

detS = 1, which diagonalizes the matrix of the operator 
5"{(f> c } and we set 



(26) 



i= 1 



where we have isolated r variables of integration 
(denoted by the tilde) which correspond to zero eigen- 
values of the operator 5"{(|) c } and do not in fact appear 
in the exponential function (25). In order to ensure cor- 
rect integration over zero-order modes, the following 
partition of unit is introduced below in the integrand (25) 



1 



= npw,-/,{<p}), 



(27) 



where \ are collective variables. An example of such a 
variable is the instanton center x defined as 



u/ d x|(p(x)| 4 (x - x ) = 0, 



(28) 



- j'tf'xoS 



(29) 



for which integration of the type (27) has the form 

( jA|q>(*)|V 
JV*x|(p(x)| 4 

By introducing collective variables (which can also be 
the instanton "orientation," its radius, and so on [5]), 
we can confine ourselves to homogeneous functions 
/{cp} [compare with (29)] where the degree of homoge- 
neity can be considered to be zero without limiting the 
generality: if/{|l(p} = \l p f {{<$>}, the substitution \ — ► 
\\PXi eliminates the factor \\P from (27). We linearize the 
arguments of the 8-functions in (27) near the saddle- 
point configuration 



1 



= f]pW,-/,{9 c }-(/;{9 c },5cp)) 



> = 1 



(30) 



= Y[\ d W(h-fi{<?c}-Jgc(fi{<?c},S<?)), 



i = 1 



and select the instanton such that X, -fi{§ c ) = (in (28) 
this corresponds to a choice of solution symmetric rel- 
ative to the point x = 0); then § c becomes a function of 
Xj. Substituting (26) and (30) into (25) and eliminating 
the 8-functions by integrating over Sep, , we have 

T -N -N+\-(M+r)l2, .rr,r, 1n 

h-\ = e g c det[f{$ c }] P 



x\^D'<p\Yl d ^c\^--M M \K) (3D 



x exp 



i[(8cp, S"{(U8cp)-W 2 ]\ 



where/' {(|) c } is an operator whose matrix consists of the 
columns /) {(|) c } and [...]/> is its projection onto the sub- 
space of the zero-order modes. The integral over D'cp 
and dt is determined by the determinant of the quadratic 
form in brackets in the exponential function (31) given 
by (-A0det[5"{(|) c }] F , although caution must be exer- 
cised when reducing this to a sum of squares [13]; the 
subscript F indicates a projection on a subspace com- 
plementary to the subspace of the zero-order modes. 
Performing elementary transformations in (31), we 
obtain a result having the form (3) where 



a - 



1 



s{<\> c y 



b = 



M+r 



c = (S{U) 



_ (M+r)/2 (27l) ( " r - r - 2)/2 det[/{(^ e }] /) 



V-det[5"{^ c }] P , 



(32) 



i = 1 



and M is the number of variables of integration con- 
tained within Dcp [this disappears from the answer on 
going over to a ratio of integrals of the type (1)]. 

Similar structural calculations can be made for the 
asymptotic form of the coefficients A K . Making the sub- 
stitution in (19) 



c? 



N 



-g 



(33) 



and isolating the dependence on JV corresponding to 
the asymptotic form (3), we have 



/*_! = (S{cM) exp(-jV + MnA0 



(34) 



where 



(35) 



x exp [ - N /A 9 ? + N -Nlng 



Setting e = 1/N and expanding F(e) as a series 

F(e) = Ao + Aie + A 2 e 2 + ... + A K e K + (36) 
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we have by analogy with (2) 

6?e (M + r-l)/2f dg 



~. f d€ (M+r-l)/2 f dg f _ (1) 



x " pi -;[i5i- 1+1 " s| - M " e 



(M) 



(37) 



The coefficients A a: are simply related to the unknown 
coefficients A K but differ from them (see below). 

For large K in (37) we can use the saddle-point 
method, for which the conditions have the form 



e„ = 



lng c 
K 



(38) 



and the function in the integrand for the saddle -point 
configuration is determined by the factor 



exp(-^-^lne ? ) 



{\ng c ) K K\. 



(39) 



Taking into account the substitution (33), the first two 
equations (38) coincide with (22) but using the solution 



4> = <t>c> 8c = 1 



(40) 



If § c (x) is a localized solution of the equation 5"{(j>} = 0, 
we know that there also exists a solution \|/ c (x) corre- 
sponding to two infinitely distant instantons (j) c (x) for 
which 5{\|/ c } = 25{(|) c }; in general a solution \|/ c can 
exist such that 5{(|) c } < 5{\|/ c } < 25{(|) c } which yields 
(43). Since the contribution (43) is larger than (41), in 
any real field theory the second of these mechanisms is 
the principal one; the first mechanism is only important 
in various nondegenerate cases such as zero-dimen- 
sional theory [see (8)] when the solution of the instan- 
ton equation is unique. 

Expanding the expression in the exponential func- 
tion (37) near the second instanton as far as quadratic 
terms in 8cp, Sg, and 8e and making the substitutions 
5g = igj, 8e = i€ c T, we have 



-K + (r+ l)/2, „ r -.-.-Mil -K 

c e c (SOJ) e 



r dt r dx r„ m (m) 



(44) 



x exp -- 



(Sep, S"{\v c }5y)---KT 2 



leads to a singularity not to a saddle point because of 
(39). Thus, other solutions of the system of the first two 
equations (38) must be sought for which two possibili- 
ties exist. 

1. Using other branches of the logarithm. In 

accordance with (38), e c is determined by the logarithm 
of g c and thus the substitution g c with integer g c — ► 
g c exp(2nis) is not an identity transformation: in this 
case we have lng c — - lng c + 2nis. For (]) = (j) c , g c = 
exp(2ra's) we have e c = 2%is/K and the contribution to 
the asymptotic form of A K is determined by the saddle 
points with s = ±1: 



The number of zero-order modes r' for the second 
instanton generally differs from r; these are isolated as 
before by introducing a partition of unit of the type 
(27), giving a dependence on K having the form 



A K ~{\n gc y K r[K+ r — r 



(45) 



In order to find the relationship between Ak and A K , we 
substitute (36) with e = 1/Af into (34), make the substi- 
tution N — - N+ 1, and use the series reexpansion rule 
given in the Appendix. As a result, we obtain 



A K ~{2%y K K\. 



(41) 



This is exactly the same mechanism as that used to cal- 
culate the corrections to the Stirling formula: the g 
dependence of the function in the integrand of (19) is 
similar to the x dependence in (1 1). 

2. Using other instantons. Let us assume that \|/ c is 
a solution of the equation S'{(])} = which differs from 
(|) c ; then on account of (38) and (39) we have the contri- 
bution to the asymptotic form 



A K = 



A K 



2%cg c 



(46) 



where c is a coefficient appearing in the Lipatov asymp- 
totic form (3) and determined by formula (32). Taking 
into account (46) and performing trivial transforma- 
tions in (44), we obtain 



A K = cA In 



r \K+ 



(47) 



A K ~K\ 



In 



S{y c y 

S{<\> c }. 



-K 



(42) 



which is larger the smaller 5{\|/ c }. The principal contri- 
bution comes from the second instanton (see end of 
Section 2) and has the lower estimate 



(ln2) K\. 



(43) 



where 



x ■ 



(27 l) ^- r '- 4)/2 det[ / '{ ¥c }] /) 



cVdet[S"{y c }], 



(r-r')/2 



(48) 
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xjnAY^^.-.vru,). 

i = 1 

In these structural calculations, we used the form of 
the functional integral (1) and the homogeneity rela- 
tions (18) typical of cp 4 theory; thus, the parameter M 
appearing in the result (32) for b was determined by the 
number of cof actors in the preexponential function (1). 
In other field theories, several fields of various types 
generally occur and the homogeneity relations differ 
from (18); nevertheless, for a wide range of problems 
the result for b has the previous form (32) but the 
parameter M has a different meaning. However, the 
parameter M does not appear in the asymptotic formula 
(47), indicating that its validity is not confined to cp 4 
theory: this is confirmed by the reasoning put forward 
in the following section. 



4. HEURISTIC DERIVATION 
OF FORMULA (47) 

In general, factorial series have an asymptotic form 
with complex parameters [see (8) and (10)] and the 
expansion coefficients I N are determined by the real 
part of some complex expression. We shall specify this 
for large N, 



I N = Re{ca F{N + b){l + A N ) 
+ ca N F(N + b) + ...}, 



(49) 



taking into account the Lipatov asymptotic form 
ca N F(N + b), the unknown power corrections to it 
denoted by A N , and the contribution of the next instan- 

ton ca N F(N + b); the corrections to the latter and the 
contribution of higher order instantons are shown by 
the ellipses. Removing the Lipatov asymptotic form 
from the brackets, we have 



x 



I N = Re\ca F(N + b) 



l+A » + cU N + 



(50) 



series (51) is much larger than the imaginary one and 
should originate from contributions which are higher in 
the hierarchy than the last term in (50); only A N can ful- 
fill this role. Combining the second and third terms in 
brackets (50), we obtain 3 



h = Re 



ca N T{N + b)\ 1 + const • £r(£ + b - b) 

L K 

(52) 



[ln( a /a) + /0] A 

A — r 

N K 



The singularity on the left-hand side of (51) is associated 
with high-order terms of the series [see the discussion of 
formula (4.10) in [5]] and the form of the general term 
given in (52) is in fact only valid for large K. Bearing in 
mind that the parameters a and b of the instanton contri- 
bution have the form (32) for a wide range of field the- 
ories, we return to the result (47) where, however, the 
coefficient Cj no longer has the specific form (48). 

We shall explain the meaning of these manipula- 
tions. As we know, the expansion of the function/(e) as 
a power series in e has a radius of convergence equal to 
the distance between the origin and the nearest singular 
point /(e) on the complex plane. For a factorial series 
the radius of convergence is zero and a singularity 
should be found for e = 0. Characteristic singularities 
generating factorial series have the form of branch cuts 
at which the discontinuity decays exponentially for 
|e| — ► [14] [see (51)]. It is deduced from the qualita- 
tive pattern established above that (a) A N has the form 
of a factorial series in l/N; (b) the divergence of this 
series is determined by the second instanton; (c) the 
contribution of the latter in (50) contains a characteris- 
tic singularity generating these series. From this it is 
logical to conclude that the second and third terms in 
brackets (50) form a single entity, being related to the 
real and imaginary parts of the same analytic function. 

This reasoning is merely based on the fact that the 
instanton contribution to the asymptotic form has the 
functional form (3). Thus, the result (47) is universal: it 
is not related to the specific field theory nor to the form 
of the quantity being studied (for example, single-par- 
ticle or two-particle Green's function). 



Assuming e = l/N, we can see that the last term has an 
intrinsic singularity for e = which may be attributed 
to the imaginary part of some factorial series [5] 



-P -a/e 



-P 



= ^Im£ W+P) i + /0 



(51) 



whose substitution into (50) leads to an expansion in 
reciprocal powers of N. It is natural to assume that the 
expression in brackets (50) is an analytic function 
whose imaginary and real parts can only appear in a 
strictly determined combination. The real part of the 



5. QUANTITATIVE RESULTS 

We shall apply these results to the problem of an 
anharmonic oscillator [9]. This can be reduced to one- 
dimensional cp 4 theory [15] in which the instantons can 
easily be investigated (in particular, by using a mechan- 
ical analogy [16]); the localized solution of the instan- 
ton equation is unique and all other solutions are 
exhausted by multi-instanton configurations containing 



3 The term A N also contains similar contributions from higher 
instantons which are small compared with those contained in 
(52). 
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several infinitely distant instantons. Thus, as \|/ c in (47) 
we need to take the two-instanton solution for which 
5{\|/ c } = 25{(|) c }, r' = r + 1 (an extraneous zero mode 
appears corresponding to the motion of two instantons 
relative to each other). Consequently, for an anhar- 
monic oscillator we have 



1 



(53) 



The dependence (53) can be compared with the results 
of Bender and Wu (5) using c 1 as the fitting parameter; 
results are plotted in Fig. 2 for c x = -1.4. 

Higher order instantons have been little studied in 
multidimensional cp 4 theory. An exception is the four- 
dimensional case for which an infinite series of instan- 
ton solutions was obtained analytically by Ushveridze 
[17]. The second instanton in this series, following the 
Lipatov one (£{<|> c } = 16?l 2 /3) has the action 5{\|/ c } = 
9ft 2 which gives the result for the asymptotic form of A K 



A K - c ; Inyg 



(54) 



(we assumed that r' = r + 3 because in view of the 
absence of spherical symmetry for the second instan- 
ton, three zero modes are added corresponding to its 
rotations in four-dimensional space). Unfortunately, 
there is no evidence that the Ushveridze series exhausts 
all the solutions; thus, the result (54) should be under- 
stood as a preliminary or lower estimate. 

A method of determining the complete series of 
higher order instantons numerically was proposed in 
[18]. It would be desirable to use this method to check 
the result (54) and to find the second instantons in all 
existing field theories. 

6. CONCLUSIONS 

Expression (4) can be used to interpolate the coeffi- 
cient function, by truncating the series at a finite num- 
ber of terms and selecting the parameters A K to ensure 
agreement with the lowest orders of perturbation theory 
known from direct diagrammatic calculations. This 
procedure is highly accurate and can reliably estimate 
the error but nevertheless is unsatisfactory in many 
respects. This is because when diverging series are 
summed, the analytic properties of the coefficient func- 
tion [19] are significant and these are reproduced quite 
incorrectly in this procedure: the coefficient function is 
assigned a multiple pole at N = but the intrinsic sin- 
gularity is lost at N = °o because of the factorial diver- 
gence of the series in (4) [see formula (51)]. Qualitative 
allowance for the functional form of the asymptotic A K 

in the form c x a x T{K + b { ) enables us to select basis 

functions exhibiting correct behavior for N — - °o 
which should give a positive effect even when the num- 
ber of fitting parameters is constant. Quantitative calcu- 




Fig. 2. Comparison between the asymptotic formula (53) for 
Cj = -1 .4 and the coefficients A K determined numerically in 
[9]. The value of |Aj | is given in [9] with a single significant 
digit (1 x 10 8 ) and the error corresponding to the range 
(0.5-1.5) X 10 8 is indicated in the figure. 

lations of the asymptotic form can be used to determine 
three parameters a x ,b x , and c 1 characterizing the coef- 
ficient function which is equivalent in efficiency to cal- 
culating the next three orders of perturbation theory. 
The calculations of a x and b x do not require functional 
integrals and can be reduced to solving nonlinear dif- 
ferential equations: the calculations of c { are of approx- 
imately the same complexity as the calculations of the 
leading Lipatov asymptotic form. This is incomparably 
easier than calculating the successive terms of a pertur- 
bation theory series where progression to the next order 
takes, on average, ten years 
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APPENDIX 



Derivation of Formula (8) 
Let us assume that two expansions exist: 

N N 2 N K 



(A.l) 



F N = 1 + tA- + 



N+P (N + pf 



+ ... + 



(N + pf 



+ ... . 



(A.2) 
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If the second series is factorial, 

B K = Re{ca K F(K + b)}, K — °°, (A3) 

it is easy to show by direct reexpansion that 

A K = R&e~ pla ca K T(K + b), K — °°. (A.4) 

Making the substitution N — - N + P - 1 in (9) and 
using (A.3) and (A.4), we obtain the result 

r(jV + p) = N^~ 1 j2nNe- N N N 

{ N (2ni) J 

and, applying it to relation (7) and using factorial series 
algebra [5], derive formula (8). 
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